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Abstract. We report here the discovery of two planet candidates as a result of our planet-search programme biased in favour 
of high-metallicity stars, using the ELODIE spectrograph at the Observatoire de Haute Provence. One of them has a minimum 
mass m2 sin; = 2.5 Mj^p and is orbiting the metal-rich star HD 43691 with period P = 40 days and eccentricity e = 0.14. 
The other planet has a minimum mass m2 sin; = 5.6 Mj„p and orbits the slightly metal-rich star HD 132406 with period 
P = 974 days and eccentricity e = 0.34. Both stars were followed up with additional observations using the new SOPHIE 
spectrograph that replaces the ELODIE instrument, allowing an improved orbital solution for the systems. 
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1. Introduction 

After the first publications suggesting the metal-rich nature 
of sta rs hosting giant planets in close orbits (lGonzaleall997 , 
Il998h . a number of works concerning the theme came out in 
the last few years jSantos et al.ll200li l2004 iFischer & Valenti 



20051: lGonzaleal2006h . With the increasing number of known 
planets, statistical studies were able to verify that the frequency 
of stars hosting a planetary companion is highly correlated with 
the stellar metallicity. The results show that the probability of 
finding a close-in giant planet is about 25-30% for the most 
metal-rich stars ([Fe/H] > 0.3) and only 3% for stars with so- 
lar metallicity. The extrasolar planet search biased in favour 
of high-metallicity stars can thus more quickly conduct to the 
discovery of planets in short-period orbits, the so-called hot 
Jupiters (f < 10 days), increasing the chances of finding plane- 
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* Based on radial velocities collected with the ELODIE spectro- 
graph mounted on the 193-cm telescope at the Observatoire de Haute 
Provence, France. Additional observations were made using the new 
SOPHIE spectrograph (run 06B.PNRCONS) that replaces ELODIE. 



tary transits. The identification of planets transiting bright stars 
provides a powerful approach to determine fundamental con- 
straints on the mechanisms of planet formation, the physical 
properties of the exoplanet, and the geometry of the system. 

Based on these assumptions, a few programmes have 
been initiated with the approach of looking for planets orbit- 
ing high- metallicity stars. One of them is the N2K consor- 
tium (Fisc her et al. which monitors nearly 2000 main- 
sequence and subgiant stars. Another project was conducted by 
our team with the ELODIE spectrogra ph at the Observatoire de 
Haute Provence (iDa Silva et al ] l2006h . From a sample of more 
than a thousand solar-type stars, we selected the more metallic 
ones after the first measurement to monitor their radial veloci- 
ties. 

Our programme has already yielded the detection of 
four hot Jupiters, with periods between 2.2 and 6.8 days 
and minimum masses between 1.0 and 2.1 Mjup, orbit- 



ing the stars HD 118203 and HD 1491430 dDa Silva et al 



' The planet aro und HD 149143 was independently discovered by 
'Fischer et al.' ('2 006() and the o ne orbiting HD 185269 was also pub- 
lished bv iJohnson et al.l j2006h . 
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2005), and HD 185269' 



2006h. HP 189733 dBouchv et al 
jMoutou et aPbOOel) . 

For the star HD 189733, additional photometric measure- 
ments have led to the observation of a planetary transit, mak- 
ing possible the determination of some paramete rs of the com- 



panion, such as mass, radius and mean de nsity (IBouchv et al 



20051: iBakos et alj|2006l: IWinn et al.ll2007h . This system is our 
best result, a very good example of what we expect to obtain 
with our programme, and revealed to be o f particular inter- 
est for further studies. Deming et al. ( 20061) analysed the in- 
frared thermal emission during a n eclipse of HP 189 733 b us- 
ing the Spitzer Space Telescope (I Werner et al. 1 12004 ') and de- 
termined the brightn ess temperature of the planet at 16 yum. 
Knutson et al.l (l2007h . doing observations with Spitzer at 8 fim., 
were able to construct a map of the temperature distribution of 
HP 189733 b, estimating a minimu m and a maximum bright- 
ness temperature at this wavelength. iFortnev & Marie vl (l2007b 
analysed the mid infrared observations of HP 189733 and sug- 
gested a possible presence of water vapor in the atmosphere of 
the planetary companion. Observations with the Hubble Space 
Telescope have also been proposed in order to perform precise 
measurements of the size and the orbital inclination angle of 
HP 189733 b (Pont et al., in preparation). 

In this paper, we report the discovery of two new planet 
candidates, resulted from our ELOPIE planet search pro- 
gramme biased towards metal-rich stars: a 2.5 Jupiter-mass 
planet orbiting the star HP 43691 with period P = 40 days, and 
a 5.6 Jupiter-mass planet in a long-period orbit ofP- 91 A- days 
around the star HP 132406. Such results are com plemented by 



additi onal measurements made using SOPHIE (B ouchv et al 
2006h . the new spectrograph that replaces ELOPIE. 



The radial velocity observations that have led to these re- 
sults are described in Sect. |2] The observed and derived pa- 
rameters of the star HP 43691 together with the orbital solu- 
tion adopted are presented in Sect. [3] The same are presented 
in Sect. |4]for the star HP 132406. In Sect. |5] we discuss the 
present and future status of the observational programme. 

2. Observations 

The stars HP 43691 and HP 132406 are both targets in our 
"ELOPIE m etallicity-b iased search for transiting hot Jupiters" 
survey (Pa Silva et al.ll2006.) . conducted from March 2004 un 



til Au gust 2006 with the ELOPIE spectrograph dBaranne et al 
11996 ') on the 193-cm telescope at the Observatoire de Haute 
Provence (France). In this programme we essentially searched 
for Jupiter-like planets orbiting metal-rich stars, assuming that 
such stars are more likely to host giant planets. 

After obtaining the first spectrum of HP 43691 and 
HP 132406, we have verified the high metallicity of these 
stars from a calibration of the surfa ce of the E LOPIE cross- 
correlation functions (Santos et al.. ,20021: iNaef ,2003) • After 
three measurements, we could clearly see in both stars a signif- 
icant radial velocity variation. We therefore conducted follow- 
up observations with ELOPIE, and we obtained 22 spectra of 
HP 4369 1 from November 2004 (JP = 2 453 333) to May 2006 
(JP = 2 453 872), and 17 spectra of HP 132406 from May 2004 
(JP = 2453 152) to June 2006 (JP = 2453 900). 



Table 1. ELOPIE and SOPHIE radial velocities of HP 43691. 
All values are relative to the solar system barycentre. The un- 
certainties correspond to the photon-noise errors. 



JD - 2400 000 


RV 


Uncertainty 


[days] 


[kms-'l 


[kms-'l 


ELODIE measurements 


53333.6255 


-29.123 


0.011 


53337.6057 


-29.098 


0.012 


53398.4118 


-29.015 


0.012 


53690.6694 


-28.945 


0.013 


53692.6430 


-28.962 


0.013 


53693.6240 


-28.983 


0.014 


53714.5653 


-28.912 


0.014 


53715.5590 


-28.903 


0.015 


53718.5549 


-28.904 


0.013 


53719.5681 


-28.885 


0.010 


53720.5323 


-28.853 


0.017 


53721.5123 


-28.892 


0.018 


53722.5237 


-28.907 


0.015 


53728.4315 


-28.969 


0.018 


53749.5333 


-28.989 


0.014 


53750.4982 


-28.958 


0.014 


53756.4444 


-28.840 


0.018 


53808.2808 


-29.091 


0.012 


53809.2850 


-29.087 


0.009 


53839.3031 


-28.960 


0.019 


53870.3405 


-28.881 


0.017 


53872.3456 


-28.923 


0.020 


SOPHIE measurements 


54044.6274 


-28.980 


0.003 


54051.6411 


-28.830 


0.004 


54053.5968 


-28.854 


0.003 


54078.6009 


-29.059 


0.004 


54079.4859 


-29.040 


0.004 


54080.4572 


-29.018 


0.004 


54081.4440 


-28.993 


0.003 


54087.4744 


-28.858 


0.004 


54088.5858 


-28.872 


0.004 


54089.6137 


-28.846 


0.004 


54142.4798 


-29.064 


0.004 


54148.4607 


-29.089 


0.004 


54151.4010 


-29.065 


0.004 


54155.4284 


-28.981 


0.004 



The ELOPIE instrument was decommissioned in August 

2006 and replaced by the SOPHIE spectrograph. Additional 
measurements were then obtained using this new instrument: 
14 spectra of HP 4369 1 from November 2006 (JP = 2 454 044) 
to February 2007 (JP = 2454155), and 4 spectra of 
HP 132406 from Pecember 2006 (JP = 2454080) to May 

2007 (JP = 2454230). 

With ELOPIE, the average signal-to-noise ratio (S/N) cal- 
culated from the spectra at /15500 A is ~40 for both stars, with a 
typical exposure time of 20 min. On the other hand, the gain in 
efficiency of SOPHIE compared to ELOPIE is more than one 
magnitude in the High-Resolution mode (used for high preci- 
sion radial-velocity measurements). Typical S/N obtained with 
SOPHIE are thus twice larger than those of ELOPIE for expo- 
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Table 2. ELODIE and SOPHIE radial velocities of HD 132406. 
All values are relative to the solar system barycentre. The un- 
certainties correspond to the photon-noise errors. 



JD - 2400000 


RV 


Uncertainty 


[days] 


[km s-'] 


[km s-'] 


ELODIE measurements 


53152.4773 


-37.821 


0.010 


53154.4825 


-37.837 


0.008 


53218.3605 


-37.858 


0.010 


53520.4238 


-37.928 


0.013 


53536.4140 


-37.875 


0.010 


53576.3681 


-37.859 


0.012 


53596.3805 


-37.818 


0.014 


53807.6666 


-37.727 


0.025 


53808.6537 


-37.755 


0.011 


53809.6643 


-37.742 


0.01 1 


53869.5117 


-37.779 


0.011 


53870.4406 


-37.771 


0.007 


53873.4386 


-37.768 


0.009 


53895.4251 


-37.757 


0.009 


53896.4387 


-37.745 


0.012 


53899.4286 


-37.743 


0.019 


53900.4376 


-37.770 


0.013 


SOPHIE measurements 


54080.7252 


-37.718 


0.003 


54173.6848 


-37.752 


0.004 


54187.6332 


-37.770 


0.004 


54230.5803 


-37.790 


0.004 



sure times 2-3 times smaller Table [T] lists the radial velocities 
of HD 43691 a nd Table[2|lists th o se of HD 132406. 

Following IZucker & Mazehl (I2OO1I) . we tried to look for 
the astrometric signatures of the two orbits in Hipparcos 
Intermediate Astrometric Data (IAD). HD 132406, whose best- 
fit Keplerian period was close to the Hipparcos mission dura- 
tion, seemed especially suitable for this kind of analysis. We 
found no evidence of astrometric signatures. Furthermore, the 
mass upper limits that the IAD produce are in the stellar regime 
and therefore do not provide any useful constraint. 

In order to derive some of the fundamental stellar parame- 
ters, like effective temperature, surface gravity and metallicity, 
using accurate spectroscopic analysis, we obtained a high S/N 
spectrum (-130 at ^5500 A) of HD 43691 with the SOPHIE 
spectrograph. 

3. A planetary companion to HD 43691 

3.1. Stellar characteristics of HD 43691 

HD43691 (HIP 30057) is listed in the Hipparcos catalogue 
(lESAll997h as a GO star in the northern hemisphere with visual 
magnitude V - 8.03, color index B - V - 0.596 and parallax 
n - 10.73 + 1.16 mas (a distance of 93 pc from the Sun ). The 
bolometric correction is BC = -0.034, derived from iFlowei 
(Il996h . Using the Hipparcos parameters we derived an abso- 
lute magnitude My = 3.18, which represents a high luminosity 
for a GO star. This suggests that HD 43691 is slightly evolved 



Table 3. Observed and estimated parameters of HD 43691 
and HD 132406. Some of the stellar parameters of HD4369I 
were obtained from spectroscopic analysis while those of 
HD 132406 come from calibrations of the ELODIE CCF 





HD 43691 


HD 132406 




Snpptral Tvnp 

vj IJL^V-' LI LlL L y 


GO IV 


GO V 




y 


8.03 


8.45 




B-V 


0.596 


0.65 




n 


10.73 ± 1.16 


14.09 ± 0.77 


[mas] 


Mv 


3.18 


4.19 




BC 


-0.034 


-0.062 




Tcff 


6200 ± 40 


5885 ± 50 


[K] 




1.38 ±0.05 


1.09 ±0.05 


Mo 


age 


2.0 - 3.6 


64 ± 0.8 


Gyr 


logg 


4.28 ±0.13 






[Fe/H] 


0.28 ± 0.05 


0.18 ±0.05 




vsin i 


4.7 


1.7 


[kms-'] 



towards the subgiant branch. iNordstrom et al. (2004) found a 
difference of 1.19 mag from the ZAJVIS, indicating the degree 
of evolution of this star 

Applying the spectroscopic analysis described in 



Santos et al 



12004) to the high S/N spectrum of HD 43691 
we obtained: T^g = 6200 + 40 K, log g = 4.28 + 0.13 and 
[Fe/H] = 0.28 + 0.05. From the calibrations of the ELODIE 
CCF (Santos et al. 2002; Naef 2003), we estimated a slightly 
smaller but compatible value for the metallicity ([Fe/H] - 0.22 
± 0.05), and a projected rotation velocity vsin / - 4.7 km s"' . 

With these stellar parameters, we estimated the mass and 
age of HD 4369 1 using the Geneva m odels of stellar evolu- 
tion computed by S chaerer et al. ( 1993 ). We found a mass of 
M* = 1.38 ± 0.05 Mq and an age between 2.0 and 3.6 Gyr, 
which are in agreeme nt with the determinations done by 
Nordstrom et al.1 (|2004 : mass = 1.38 + 0.08 Mq and age 
2.6 ± 0.5 Gyr. These values are compatible with the star be- 
ing slightly e volved, especially ta king into account the stel- 
lar metallicity iMowlaviet al.ll998 ). The observed and derived 
stellar parameters of HD 4369 1 are shown in Table[3j 

3.2. Orbital solution for HD 43691 b 

The best Keplerian orbital solution fitted to the radial velocities 
of HD 43691, using both ELODIE and SOPHIE observations, 
provides a orbit with period P = 36.96 + 0.02 days and ec- 
centricity e = 0.14 + 0.02. With the estimated value for the 
primary mass of 1 .38 we obtained a minimum mass m2sin / 
- 2.49 Mjup and a separation of 0.24 AU for the planetary com- 
panion. The solution includes a velocity zero-point of the two 
datasets as a free parameter, and the difference between them is 
As-E = 23 ± 4 ms"'. For this solution, we estimated a false 
alarm probability of 1 .3 x 10"^ using the approach described in 
Horne & Baliunasl(ll986l) . 

In the top panel of Fig. [T] we plot the radial velocities of 
HD4369I and the Keplerian fit adopted using the two sets of 
measurements. The middle panel shows the residuals around 
the solution. The weighted rms around the solution is cte = 
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Table 4. Orbital elements for the best Keplerian solution of 
HD 43691 and HD 132406 as well as the inferred planetary pa- 
rameters. 
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Fig.l. Top: ELODIE and SOPHIE radial velocities of 
HD 43691 plotted together with the best Keplerian solution that 
fits the combined measurements. Middle: Residuals around the 
solution, with cte = 17.5 ms-' for ELODIE, erg - 9.0 m s ' for 
SOPHIE, and ctes - 10.0 ms ' for the combined data points. 
Bottom: Phase-folded radial velocities with the best Keplerian 
solution. Error bars represent the photon-noise errors. 



17.5 ms ' for ELODIE, os = 9.0 ms ' for SOPHIE, and 
Ces = 10.0 ms ' for the whole dataset. The bottom panel 
shows the phase-folded radial velocities. Table H] lists the 
adopted orbital elements, together with the inferred planetary 
parameters. 

3.3. Low chromospheric activity for HD 43691 

The radial velocity variations observed for a star can also be the 
result of physical events in the stellar atmosphere rather than 
the presence of an orbital companion. For example, spots on the 
surface of an active star can change the observed spectral-line 
profiles and induce periodic variations of the measured radial 
velocities. By analysing the line-bisector orientations one can 
distinguis h which of these situations is the real origin of the 
variations (lOueloz et al..,2001.) . 





HD 43691 


HD 132406 




p 


36.96 + 0.02 


974 ± 39 


[days] 


T 


54046.6 + 0.5 


53474 + 44 


[JD - 2400000] 


e 


0.14 + 0.02 


0.34 ± 0.09 




V 


-29.000 ± 0.003 


-37.840 ± 0.008 


[kms-'] 


a> 


290 ± 5 


214 ± 19 


[deg] 


K 


124 ± 2 


115 ± 26 


r —It 

[ms '] 


^meas 


22 (E) + 14 (S) 


17 E + 4 S 




As_E 


23 ± 4 


93 ± 17 


[m s-'l 




17.5 


12.1 


[m s-'] 


o-s 


9.0 


4.1 


[ms-'] 


0"ES 


10.0 


7.5 


[ms-i] 


flisin / 


4.17 


9.73 


[10-" AU] 


f{m) 


7.06 


1.30 


[10-" Mo] 


ra2sin i 


2.49 


5.61 


[Mj„p] 


a 


0.24 


1.98 


[AU] 



The analysis of the line-bisector orientations, or bisector 
inverse slope (BIS value), of HD 43691 shows that there is no 
correlation between the BIS values and the radial velocities de- 
rived (Fig. |2] top panel). Thus the observed variations in ra- 
dial velocity are not induced b y stellar activity an d rotation 



(as is the case of HD 166435 in lOueloz et al.ll2001h . The ob- 
served behaviour of the line bisectors also indicates that the 
radial velocity variations are not resulting from contamination 
by the light o f a late-type binary companion (see e.g. the case of 
HD 41004 in lSantos et al.ll2002l) . Such variations are thus most 
probably due to the presence of a planetary companion orbiting 
HD 43691. 

The chromospheric activity level can also be verified by 
means of the reemission in the core of Can absorption lines 
(e.g. /13968.5 A). By observing the respective spectral region 
in the high S/N spectrum of HD 43691 obtained with SOPHIE 
(Fig. |2] bottom panel), we can note that this star is not active. 
Since this line is located in the blue part of the spectral domain, 
where the flux is appreciably lower, this kind of inspection re- 
quires a high S/N spectrum, which is much better achieved by 
the higher efficiency of SOPHIE. 

4. A long-period planet orbiting HD 132406 

4.1. Stellar characteristics of HD 132406 

HD 132406 (HIP 73 146) is listed in the Hipparcos catalogue as 
a GO star with visual magnitude V = 8.45, color index B -V - 
0.65 and parallax n - 14.09 + 0.77 mas (71 pc distant from the 
Sun). These parameters set a value of My = 4.19 for the abso- 
lute magnitude. The bolometric correction is BC = -0.062. 

The metallicity and the projected rotation velocity of this 
star are respectively [Fe/H] - 0.18 + 0.05 and vsin/ =1.7 
kms-', estimated from the calibrations of the ELODIE cross- 
correlation functions. The efifective temperature derived is Tes 
- 5885 + 50 K and comes from the cafibrations of T^s as a 
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-29 
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-28. ( 
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Fig. 2. Top: Comparison between bisector inverse slope (BIS) 
and radial velocities of HD 43691 showing no correlation be- 
tween them for both ELODIE and SOPHIE measurements. 
Bottom: /13968.5 A Ca ii absorption line region of the high S/N 
spectrum obtained for HD 43691. No clear emission feature is 
observed in the center of this line, indicating a low activity 
level. 



function ofB-V and [Fe/H] fromlSantos et all (|2004 . The de- 
rived mass and age are M* = 1 .09 + 0.05 Mq and 6.4 + 0.8 Gyr, 
from the Geneva models of stellar evolution. These parameters 
are listed in Table [3] 

4.2. Orbital solution for HD 132406 b 

A Keplerian solution fitted to the radial velocity measurements 
of HD 132406, from both ELODIE and SOPHIE observations, 
results in an orbit with period P = 974 + 39 days and eccen- 
tricity e - 0.34 + 0.09. The velocity zero-point of the two 
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53500 54000 
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Fig. 3. Top: ELODIE and SOPHIE radial velocities of 
HD 132406 plotted together with the adopted Keplerian so- 
lution that better fits the combined measurements. Bottom: 
Residuals around the solution, with cte = 12.1 ms"' for 
ELODIE, crs = 4.1 ms ' for SOPHIE, and ctes = 7.5 ms ' 
for the combined data points. Error bars represent the photon- 
noise errors. 



datasets is a free parameter, and the diff'erence between them 
isAs-E = 93 + 17 ms-'. 

The top panel of Fig. |3] shows the radial velocities of this 
star together with the best Keplerian solution. The bottom panel 
of the same figure shows the residuals around the adopted so- 
lution, for which the weighted rms is cte - 12.1 ms ' for 
ELODIE, crs =4.1 ms ' for SOPHIE, and ergs = 7.5 ms"' 
for the combined set of measurements. HD 132406 is slightly 
fainter than HD 43691, but has smaller photon-noise errors, 
which is probably due to broader line profiles of HD 43691. 
Table |4] lists the orbital elements and the planetary parameters 
of the HD 132406 system, which was obtained with the best 
Keplerian fit. 

As in the case of the star HD 43691, the analysis of the 
bisector inverse slope of HD 132406 shows no correlation be- 
tween the BIS values and the observed radial velocities (Fig.|4]l. 
In addition, no chromospheric reemission is observed in the 
core of the Ca ii absorption line at /13968.5 A. 



5. Discussion and Conclusions 

In this paper we have announced the discovery of two new 
planet candidates as results from our ELODIE search pro- 
gramme biased towards metal-rich stars. In this programme, a 
total of six planets were discovered so far, out of which four are 
hot Jupiters (f < 10 days) and two are the intermediate-period 
planets just presented. Five of the host stars have metallicity 
greater than 0. 1 dex, while the one with a transiting very hot 
Jupiter (HD 189733) is a solar-metallicity star. 
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ELODIK o SOPHIE 



s 




Radial velocity [km s '] 

Fig. 4. Comparison between bisector inverse slope (BIS) and 
radial velocities of HD 132406 showing no correlation between 
them. Both ELODIE and SOPHIE measurements are plotted. 

So far, we have observed at least once almost 82% of the 
1061 sample stars, and for each one we estimated a value for 
metallicity. We thus verified that about 26% of the observed 
stars have [Fe/H] > 0.1 dex (about 15% for 0.1 < [Fe/H] < 0.2, 
8% for 0.2 < [Fe/H] < 0.3 and 3% for 0.3 < [Fe/H] < 0.4). 
Furthermore, according to the pe rcentage of s tars with planets 
per metallicitv bin determined bv lSantos et al. (2004), 10-30% 
of the stars with [Fe/H] > 0.1 dex are likely to host a giant 
planet (about 9, 24 and 28% respectively for the same three 
ranges of metallicity mentioned above). Applying these per- 
centages to the 867 observed stars, we finally find that the 
number of giant planets we expect to discover in each of those 
metallicity ranges is respectively 12, 17 and 7, a total of 36 
planets, from which roughly 25% (9 planets) are predicted to 
be hot Jupiters. 

Although most of our target stars were already observed, 
only 75% of the metal-rich stars have a minimum of three mea- 
surements. Stars with one or two spectra need more observa- 
tions before being rejected or classified as possible planet hosts. 
On the other hand, long-period planets are more difficult to de- 
tect, and stars showing long-term radial-velocity trends also 
need more observations. In any case, we have akeady found 
almost a half of the expected number of hot Jupiters among the 
metallic portion of our sample. 

The stars presented in this paper were also observed using 
the new SOPHIE spectrograph. The proposed orbital solutions, 
firstly found with ELODIE, were improved with the new obser- 
vations. With ELODIE decommissioned, this new instrument 
will also continue monitoring other high-metallicity stars, es- 
pecially the most promising cases. 
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